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^ ■ Abstract 



This study uses a combination of stochastic optimization, statistical mechanical theory, and 



, molecular simulation to test the extent to which the long-time dynamics of a single tracer parti- 
es 

Q ' cle can be enhanced by rationally modifying its interactions-and hence static correlations~with 
the other particles of a dense fluid. Specifically, a simulated annealing strategy is introduced 
that, when coupled with test-particle calculations from an accurate density functional theory, 

o ; 

T— I , finds interactions that maximize either the tracer's partial molar excess entropy or a related 

in 

CN ' pair-correlation measure (i.e., two quantities known to correlate with tracer diffusivity in other 
contexts). The optimized interactions have soft, Yukawa-like repulsions, which extend beyond 
I , the hard-sphere interaction and disrupt the coordination-shell cage structure surrounding the 
tracer. Molecular and Brownian dynamics simulations find that tracers with these additional 

X' 

^ • soft repulsions can diffuse more than three times faster than bare hard spheres in a moderately 
supercooled fluid, despite the fact that the former appear considerably larger than the latter by 
conventional definitions of particle size. 



^ These authors contributed equally to this work. 
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I. INTRODUCTION 



Diffusion of a tracer particle through a complex fluid is an essential process for many 
chemical, materials, and biological systems. The rate at which it occurs-characterized, e.g., 
by the corresponding long-time diffusion coefficient-has practical implications for techno- 
logical developments in drug delivery, catalysis, materials synthesis, and separations.-"- 
Experimental studies clearly show that there is an intimate connection between tracer 
diffusivity and the interparticle interactions present in a system.-"- Furthermore, ap- 
proximate microscopic approaches for predicting dynamics, such as kinetic,-i^ generalized 
Langevin,— "— and mode-couplingi^"— theories, have recently provided some important in- 
sights into this link. Nonetheless, a comprehensive theoretical framework-or even a clear 
conceptual picture-for understanding and predicting how tracer diffusivity should depend 
on the interparticle interactions is still lacking. As a result, basic questions about how to 
engineer complex fluids with desired dynamic properties remain challenging to address, 
especially for cases where nontrivial constraints on particle properties (size, chemistry, 
etc.) and thermodynamic conditions (e.g., particle concentration, composition, etc.) must 
be satisfied. 

In this paper, we present theoretical results that shed new light on basic physics relevant 
to designing single-particle dynamics of complex fluids. In particular, we explore-within 
a model system-the extent to which one can rationally modify the tracer diffusivity of a 
particle by tuning how it interacts with the other particles in the system. For colloids or 
nanoparticles suspended in solvent, this is particularly interesting because it is often pos- 
sible to systematically modify their effective interparticle interactions by varying solvent 
quality (e.g., by changing composition of a third component such as a surfactant or salt), 
by introducing a depleting agent (e.g., a nonadsorbing polymer), by chemically attaching 
or physically adsorbing molecules onto the particles, or by tuning of electric or magnetic 
interactions via an external field.—"— 

As a starting point for our analysis, we consider a single hard-sphere tracer particle in 
a dense fluid of other hard spheres. The question we address here is simply, which types 
of isotropic pairwise potentials between the tagged tracer particle and the neighboring 
particles-when added to the bare hard-sphere interaction-significantly increase the tracer's 
long-time diffusion coefficient? 
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At first glance, the idea that adding a contribution to the tracer particle's hard-sphere 
interaction could significantly increase its long-time mobility may not be intuitively ob- 
vious. The naive expectation is that adding an attractive interaction would slow down 
the dynamics of a single tracer because it would energetically stabilize the surrounding 
"cage" structure formed by its nearest neighbors. On the other hand, adding a repulsive 
interaction to the pair potential might likewise be expected to decrease the mobility of 
the tracer because it would effectively increase its size. Based on these arguments, one 
might hypothesize that the bare hard-sphere interaction is optimal for tracer diffusivity. 
Interestingly, as we discuss in detail below, adding certain types of softer, repulsive inter- 
actions to the tracer's hard core can result in a pronounced increase in its mobility. The 
key to achieving the enhanced diffusivity appears to be finding interactions that disrupt 
the coordination-shell structure surrounding the tracer without significantly increasing the 
tracer's effective size. 

How can one discover the mathematical forms of such mobility enhancing interactions in 
a reasonably efficient and systematic way? The strategy we pursue here is to first identify 
structural or thermodynamic quantities, relatable to interparticle interactions through 
equilibrium statistical mechanics, that correlate strongly with tracer diffusivity in other 
related systems. Once an appropriate static property has been identified, we employ a 
variational approach that uses either computer simulations or liquid-state theory to find 
the types of interparticle interactions that increase (perhaps even optimize) the static 
property, subject to relevant realizability constraints. Finally, we perform molecular or 
Brownian dynamics simulations to determine whether the interactions do in fact result 
in tracer-particle long-time diffusivities significantly higher than those of the reference 
hard-sphere particles. 

The specific strategy we adopt here is motivated by the observation from molecular 
simulations and experiments that the excess entropy (relative to an ideal gas in the same 
volume) s*^^, as well as its associated two-particle approximation s^'^\ can be used to 
semi- quantitatively predict a number of nontrivial effects that not only temperature and 
particle concentration^^-^but also confinement and other external fields^^'^^^-^have on 
the transport coefficients of equilibrium and moderately supercooled fiuids. Excess entropy 
is a static measure that characterizes the reduction in the number of states available to a 
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system due to interparticle correlations. Although a rigorous link between excess entropy 
and dynamics that can explain all of the aforementioned observations is still lacking, 
the connection is now well understood for low-density gases,— fluids with inverse-power- 
law pair potentials, and fluids that exhibit so-called "isomorphs" M It is also physically 
reasonable to expect that the collisional processes important for dynamics correlate with 
structural measures that track the strength of the static interparticle correlations (at least 
for conditions not too close to the glass transition). 

There are two excess-entropy based quantities that are natural to investigate when 
the focus is on the structure and dynamics of a single tracer particle: (i) the two-body 
contribution to excess entropy, si , arising from the partial radial distribution function 
between the tracer and its neighbors, and, (ii) the partial excess molar entropy of the 
tracer sf^. Both have been found in molecular simulations to correlate with the long-time 
dynamics of particles in mixtures.— 

We organize the balance of the paper as follows. In Sec. [Tll we provide information 
about the model system that we study, and we discuss how to compute si and sf^ for 
this system from statistical mechanics. We then discuss simple variational schemes for 
obtaining tracer-fluid potentials which result in significantly higher values of si and sf^ 
as compared to those associated with the bare, hard-sphere tracer. We also present details 
concerning how we compute the corresponding tracer diffusivities via molecular (and, in 
a few cases for comparison, Brownian) dynamics simulations. In Sec. IIIIl we show the 
link between tracer diffusivity and the two aforementioned static measures for hard-sphere 
tracer particles of various size. We then demonstrate that the variational calculations 
discussed above produce soft, Yukawa-like repulsive potentials that extend beyond the 
bare hard-sphere interaction and disrupt the coordination-shell cage structure surrounding 
the tracer. Particles with the additional soft potentials can have tracer diffusivities more 
than three times higher than those of the bare hard-sphere particles in the moderately 
supercooled fluid. This is striking, especially given that the former would be considered 
"larger" than the latter by conventional measures (e.g., Barker-Henderson radius, second- 
virial coefficient, partial molar volume, and solubility). Finally, we present calculations 
using a model tracer particle with hard-core plus Yukawa interactions that illustrate how 
the optimized interactions presented in this study strike a balance between being strong 



4 



enough to reduce local coordination-shell structure that impedes diffusivity, but still soft 
enough to avoid typical decreases in mobility due to larger particle size. 



II. THEORY AND METHODS 



A. Model Interactions 



Consider a single tracer (t) particle in a bath of hard-sphere-like fluid (f ) particles. We 
model the hard cores of all particles in the system via pair potentials of the following 
Weeks- Chandler- Andersen (WCA) form,— 



iWCA, 



r;d,e) 



r > 2i/24rf 



The quantities d and e characterize the usual length and energy scales of the WCA inter- 
action, respectively. The fluid particles surrounding the tracer have a nominal diameter a 
and interact solely via this potential, 

</'ff(r) = 0^^^(r;a,6) (2) 

The pair interaction between the tracer particle and the surrounding fluid particles is given 
by: 

0tf(r) = < (3) 

Here, cxtf = (crt -f- o")/2, and dt is the nominal diameter of the tracer particle. The addi- 
tional interaction, 4>Q{r), is the contribution that can be "tuned" in order to maximize s\ 
or s^^, and -as we show- increase the dynamics of the tracer particle. The procedure for 
computing (pair) is discussed in Sec. Ill Dl and IIII A[ To ensure that the tuning procedure 
does not impact the size of the exclusion region that constitutes the bare hard-core in- 
teraction, we define tq such that gtii^) > 5 x 10"'^ for r > rg (and gtii^) < 5 x 10^^ for 
r < tq), where gti{f) is the partial radial distribution function (PRDF) between the tracer 
and the other fluid particles. For the range of parameters studied in this work, this leads 
to ro = (0.97-0.98)atf, and <p^^^{ro] a^u e) = S-S/crT. 
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B. Calculation of s^^^ and s^^ 

Consider a binary fluid mixture of A^t tagged tracer particles and N other particles 
in the inflnite dilution limit of the former [i.e., vanishing tracer-particle mole fraction 
Xt = Nt/{Nt + N) — > 0]. The two-body contribution to the excess entropy of the fluid 
arising from the tracer-fluid positional correlations in this limit is given by^i^ 

(2) pkB 



2 



{gti{r) lngti{r) - gti{r) + l}dr (4) 



where p = N/V, V is the volume, and k-g is Boltzmann's constant. 

More generally, the partial molar excess entropy associated with the tracer particle is 
given by 



dS \ fdS 



(5) 



where 5, T, and P are the entropy, temperature, and pressure of the fluid, respectively. 
S^^ is the corresponding entropy of an interaction- free (i.e., ideal gas) version of the system 
with the same values of V , Nt, and A^. Using standard thermodynamic relations, one can 
recast sf^ as follows, 

if /fcB =PW+{Z- l)pv, - /3pT (6) 

Here /3 = (k^T)'^, and Z = 13 P/ p is the compressibility factor of the fluid. The quantities 
e^^, Vt, and pl^ represent the tracer particle's excess energy, partial molar volume, and 
excess chemical potential, respectively. 

For the infinitely dilute tracer particle in a fluid of hard-spheres considered here, each 
term on the right-hand side of eq. Elcan be simplifled considerably. Speciflcally, the partial 
molar excess energy of the tracer is given by 

er = P j gtfir)(j)tf{r)dv (7) 

The partial molar volume for the tracer can be calculated using the Kirkwood-Buff 
relationship:— 

Vt = KkBT + j [1 - (?tf(r)]rfr (8) 

where k is the isothermal compressibility of the pure hard-sphere fluid, which we estimate 
here using the Carnahan-Starling equation of state.^- Note the integral in eq. |8] converges 



slowly and hence its numerical value is sensitive to the truncation point. However, as 
discussed elsewhere,— its value can be accurately estimated by extrapolating the values 
obtained when integrating out to various maxima and minima of gt{{r) (in the range 
r = 6 - 10(t). 

Finally, the excess chemical potential of the tracer is equal the reversible work of in- 
serting a test tracer particle into the hard-sphere fluid at constant temperature T and 
chemical potential fj,,— 

/ir = fitMr);T,/i]-fi(/i,T) (9) 

Here, Qt is the grand-potential of the fluid with the tracer particle fixed at the origin, 
and Q is the grand-potential of pure fluid without the tracer particle. Density functional 
theories for inhomogeneous fluids, like the one discussed in Sec. Ill C[ provide an approx- 
imate functional relationship for Qt[gti{r);T, fi], while Q{fi,T) is a property of the pure 
hard-sphere fluid. In this work, the latter is again estimated from the Carnahan-Starling 
equation of state.— 

A key observation of this section is that, for an infinitely dilute tracer in a hard-sphere 
fluid, s[^^ is a functional of gti{r), and sf^ is a functional of gt{{r) and (f)tf{r). Moreover, 
a basic result of statistical mechanics^ is that gt{{r) uniquely determines (j)ti{r)-a.nd vice 
versa-for this problem. As a consequence, one can search out optimal functions for gtf{r) 
[or equivalently 0tf(^)] that maximize either si or sf^ subject to suitable realizability 
constraints. 

C. Density functional theory 

Classical density functional theory (DFT) of inhomogeneous fluids provides a conve- 
nient means for interconverting between gt{{r) and 4>ti{f)- In this work, we use a recent 
modification^ of Rosenfeld's fundamental measure theory^>^ which-by construction- 
reproduces the Carnahan-Starling equation of state for a hard-sphere fluid in the homoge- 
neous limit. For our DFT calculations, we map the fluid of hard-sphere-like WCA particles 
onto "equivalent" hard spheres of diameter cTf, an approximation that has been found to 
be adequate in other contexts.— The tracer particle, on the other hand, is treated as a 
fixed external potential of the form (ptfi^) [see eq. [3]. This "test-particle" approach allows 
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us to use the inhomogeneous one-component version of fundamental measure theory to 
compute infinite-dilution partial molar properties of the tracer particle in a homogeneous 
hard-sphere fluid.— 

The following expression provides a convenient representation for computing gti{r) 
within this framework, 

= ^exp{P[fi-Mr)]+c'^'\r,g^i{r);fi]} (10) 

where Ab is the thermal deBroglie wavelength of the fluid particles and c^^^ is their one- 
body direct correlation function in the external field of the test particle. The latter can 
be expressed as a functional derivative of excess intrinsic Helmholtz free energy J-"^^ with 
respect to gtf{r), 

Ji) _ UJ^""[^tf(r);/i] 

P Sg,,{r) ^''^ 

Explicit functional expressions for both J-"*^^ and Qt of eq. [9] from the modified fundamental 
measure theory are presented elsewhere.— 

To obtain gtf{r) for a given 0tf(^); T, and fi, we solve eq. [10] using Picard iterations.— 
Alternatively, for a given gtf{r), T, and fi (and hence p), we obtain the corresponding 
4'ti{r) by numerically inverting eq. [TOl 



D. Optimizing tracer-fluid interactions 



1. Optimizing s[^^ 



From eq. HJ it can easily be seen that a "flattened" PRDF-i.e., gti{r) = 0(r — (Xtf), where 
0(r — (Ttf) is the Heaviside step function-maximizes si subject to the hard-core constraint. 
We insert this target gti{r) into the inverse DFT protocol described in section III CI to 
determine the corresponding tracer-fluid interactions 0tf('") for various values of (Xt/cr and 
reduced fluid density pa^. We also compute the partial molar excess entropy sf of this 
particle with the methods described in section III Bl 
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2. Optimizing sf^ 

To compute tracer-fluid interactions 0tf(^) [and corresponding Qtii^)] that minimize 
—sf^/k-Q, subject to the hard-core constraint, we use a simulated anneahng global op- 
timization algorithm described in Corana et al.— . For the optimization algorithm, we 
assume that 5'tf(^) for '"o < ^ ''"a can be approximately described by the following ex- 
pression: 



Here, Am = "ra-njij^ — ro), and is chosen such that \gti{r^ — 1| < 0.02 for a single 
hard-sphere-like tracer particle of diameter at surrounded by a fluid of hard-sphere-like 
particles of diameter a with density pa^ . For r > ra, we adopt a standard asymptotic 
form for guir^^ and determine its required coefficients from a least-squares fit to the same 
hard-sphere tracer's PRDF data. 

The optimization algorithm proceeds as follows. An initial coefficient vector 
[ai, 035, 61, 65] is chosen by fitting eq. [12] to simulation data for the equilibrium 
structure surrounding the aforementioned hard-sphere-like tracer particle. The objective 
function [— s'^t' ^^^^cnJ f*^^ ^^^^ system is then computed using the method described in 
section IIIBI An initial effective temperature for the simulated annealing algorithm, k^T^^ 
(= 1 X 10~^ in this work) is also provided as input. 

A cycle in the simulated annealing scheme is then carried out as follows, (i) A "trial" 
coefficient vector is created by randomly modifying the coefficients of the "current" coef- 
ficient vector. The magnitudes of the modifications are chosen such that approximately 
half of trial coefficient vectors are accepted according to the criteria described in the next 
step, (ii) If the trial gtiix^ does not satisfy the constraints described in section Hi At then it 
is rejected and step (i) is repeated. Otherwise, it is input into the inverse DFT routine of 
section Hi CI to determine (^i,i(r\ and the corresponding partial molar excess entropy s^^triai 
is computed as described in section III CI and IIIBI respectively. If the trial guir^ results 
in — St^'^j.jai < then it is accepted and it replaces the current PRDF. If not, the 

trial PRDF is accepted with probability exp[(s5;''| trial — ■St''| current )/^B^sa]- In the present 
study, approximately 2 x 10^ trial moves consisting of steps (i) and (ii) are carried out per 




(12) 
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cycle. 

After completion of a cycle, a lower /cB^sa is chosen that is 85% the value of the previous 
simulated annealing temperature. A cycle at this lower temperature is then initiated unless 
the run has converged according to the following metric: the difference in the minimum 
value of —sf^/k-Q obtained in four successive temperature reduction cycles is less than 0.1. 

E. Dynamics simulations 

Molecular dynamics simulations were carried out in the microcanonical ensemble using 
the velo city- Ver let integration method with a time step of O.OOlcry/m/e. The average 
temperature of all runs was e/fce, which was set by periodic velocity rescaling during pre- 
equilibration simulations.^ For comparison, a small number of Brownian (overdamped 
Langevin) simulations were also carried out using the Ermak algorithm using timestep 
At = 0.002rB where tb = mDo/e and Dq = O.OOIcta/^/c.^^ A single tracer particle 
and 4000 < N < 8000 hard-sphere-like fluid particles per (periodically-replicated) cubic 
simulation cell of volume V were used; the values of these parameters in a given simulation 
were chosen to realize a specified reduced density pcr^ of the fluid. The tracer diffusivity 
was computed by fitting the long time behavior of the mean squared displacement to the 
Einstein relation for diffusion (Ar^) = 6Dt. Error bars represent 95% confidence intervals 
determined by analyzing the results of between 30-192 independent trajectories. 

III. RESULTS AND DISCUSSION 

In this section, we use methods described above to compare how the static structure 
and single-particle dynamics of an infinitely dilute hard-sphere-like tracer particle in a 
dense hard-sphere fiuid compare to the same quantities for tracers whose hard-sphere 
interactions are augmented by additional pair interactions to optimize either si or sf^. 
Except where indicated, all results for PRDFs and long-time diffusivities are obtained 
from molecular dynamics simulations. 

As a starting point. Fig. [T] shows pair potentials and PRDFs for the case of (Xt = 2cr and 
pcr'^ = 1. This is a moderately supercooled fluid; the thermodynamic freezing transition 
for hard spheres occurs at pa^ ~ 0.94. The main feature to note in Fig. [1^ is that the 
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FIG. 1: (a) Potential [</>tf(?')] and (6) PRDF [s'tfC^)] between tracer particle of diameter at = 2(t 



optimized (blue dashed), and sjr^-optimized (red dot-dashed) interactions are shown. 

optimized tracer-fluid potentials differ considerably from the hard-sphere-like form due to 
the presence of longer-range soft repulsions. What effects should these repulsions have 
on tracer-particle diffusivity? On one hand, one might expect-based on arguments from 
hydrodynamics^^^that the optimized tracers will show slower dynamics due to their larger 
effective size. On the other hand. Fig. lb clearly illustrates that the main effect of the 
optimized potentials is to disrupt the coordination-shell "cage" structure surrounding the 
optimized tracers. Based on this disruption, we expect the diffusivity to increase. 

Fig. |2K shows that the diffusivity of tracers with optimized potentials has indeed 
increased by up to a factor of three (two) for the sf- {s\ -) optimized potentials, respec- 
tively. For a hard-sphere tracer, increasing the fluid density results in more pronounced 
coordination-shell "cages" of surrounding fluid particles. The s\ and sf^ optimization pro- 
cedures work to destroy this caging structure at any density, which is why the structural 
and dynamic consequences of optimization are more pronounced in denser fluids. Figj2]D il- 
lustrates that the difference between the two types of optimized potentials, clearly reflected 
in the tracer dynamics, does not manifest in s{: . This is consistent with the observation 
that the shape of the pair correlation function alone cannot generally account for the dy- 
namics of particles in dense fluids.— However, other aspects of static structure capture 



and fluid particles of diameter a and density pa' 
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FIG. 2: (a) Ratio of diffusivity (D) for si - optimized (black square) and sf^ - optimized 
(red circle) tracer to that (Dhs) of a hard-sphere-like tracer with the same exclusion diameter 
(fjt = 2a) in a fluid of hard-sphere- like particles of diameter a and density po"^. Enhancement of 
(b) si and (c) sjr^ (over hard-sphere-like tracer values) due to optimization. 

the difference. As is shown in Fig|2l3, when comparing tracers with the same hard-core 
diameter, those with higher diffusivities have considerably higher values of sf^. 

The effect of si - and s^^-based optimizations for tracer particles with different bare 
hard-core diameters (relative to those of surrounding fluid particles) is explored in Figures 
[3] and m respectively. For an infinitely large tracer, the first coordination shell of fiuid 
particles constitutes an insignificant fraction of the tracer particle size. Any modification 
to this structure is expected to have a negligible effect on tracer diffusivity.^^ At the 
other extreme ((Tt/o" = 0), there is less structuring surrounding the point-like tracer to 
begin with, and thus the benefits of eliminating this structure through optimization are 
outweighed by the associated increase in effective particle size. Between these two limits, 
there is a maximum in the diffusivity ratio (as a function of tracer particle size) for si 
and sf^ optimized potentials. In the next section, we introduce a model "hard-sphere + 
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FIG. 3: (a) Ratio of diffusivity (D) for sj: - optimized tracer to that (Dhs) of a hard-sphere- 
like tracer with the same exclusion diameter, at- Both tracers are in a fluid of hard-sphere-like 
particles of diameter a and density po"^ = 1. (6) The corresponding differences in and sf^ 
between the si - optimized tracer and hard-sphere-like tracer. 

Yukawa" tracer to investigate the dynamic consequences of the trade-offs between effective 
particle size, softness of the interparticle potential, and fluid structure. 



A. Tracer-fluid interactions with increased and s^^: Mapping onto effective 
potentials 

The form of the soft, repulsive potential resulting from optimizing sf^ looks similar to 
a HCY {i.e., hard-core [^'^'^'^(r)] -|- Yukawa interaction, which is often used 

to approximately model particles with repulsive, screened-electrostatic forces. Separating 
the Yukawa contribution from the hard-sphere contribution allows us to tune in the soft 
repulsion in a controlled way via the parameter a, where a = 1 corresponds to a best-flt 
of the model to the s5;^-optimized potential for at = 2a and density pa^ = 1. 

0HCY(^) ^ 0^t^^(r) + e~<^~''^ (13) 

r 

The parameters for the flt are ey = 13.41e and na = 0.82. 

Fig. |5] shows HCY tracer pair potentials and PRDFs for various values of a. Moving 
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FIG. 4: (a) Ratio of difFusivity (D) for sjr^ - optimized tracer to that (Dhs) of a hard-sphere- 
like tracer with the same exclusion diameter, at- Both tracers are in a fluid of hard-sphere-like 
particles of diameter a and density pcr^ = 1. One data point (green diamond) is from Brownian 
dynamics simulations described in the text. (6) The corresponding differences in si and s^^ 
between the sf^ - optimized tracer and hard-sphere-like tracer. 

from a = to a = 1, the main effect on the PRDF is to decrease the height of the first 
peak and reduce the structuring in more distant coordination shells. Fig. [6^ shows a 
corresponding enhancement of the tracer diffusivity, which might be expected given that 
the potential and structural changes resemble those of the s^^-based optimization process. 

However, for a > 1, the model tracer- fluid potential exhibits a Yukawa repulsion that 
is stronger than optimal. In fact, as can be seen from the PRDF in Fig. |5]d, the resulting 
tracer exclusion core [separations for which gtf{r) ~ 0] becomes noticeably larger than that 
of the underlying WCA potential. Accordingly, Fig. |6] shows that the tracer diffusivity 
decreases with increasing a for a > 1, as should be expected since the main effect in this 
range is to increase the effective hard-core diameter. Note that such potentials with differ- 
ent hard-core diameters than the underlying WCA potential are avoided by construction 
in the s^^-based optimization. 
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FIG. 5: ((a) Potential [(ptiir)] and (6) PRDF [gtf{r)] for the HCY tracer particle with at = 2a 
in a fluid of hard-sphere-like particles of diameter a and density pa'^ = 1. Cases a = (green 
solid), a = 0.5 (black dashed), a = 1 (blue long-dash), and a = 2 (black dot-dashed) are shown. 
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FIG. 6: Ratio of diffusivity (D) for the HCY tracer to that (Dhs) of a hard-sphere-like tracer 
with the same exclusion diameter (at = 2a) in a fluid of hard-sphere-like particles of diameter a 
and density pa^ = 1. 

IV. CONCLUSIONS 



Using a combination of simulations and liquid-state theory, we have shown that the 
long-time diffusivity of a hard-sphere-like tracer particle with a given exclusion diameter 
in a dense hard-sphere fluid can be significantly enhanced by adding a soft repulsion to 
its interactions with neighboring particles. An important effect of this repulsion is to 
disrupt the coordination shells that otherwise surround the bare hard-sphere tracer. The 
form of the required repulsions can be determined theoretically by maximizing either the 
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tracer's partial molar excess entropy or its partial radial-distribution function contribution 
to the fluid's excess entropy. We show that the long-time diffusivities of tracers with s^^- 
and St^'-optimized interactions can be higher than those of the corresponding hard-sphere 

tracers by more than factor of two and three, respectively, for moderatately supercooled 
liquid-state points. 

In future studies, we plan to investigate the effect of the aforementioned soft repulsions 
on the position-dependent diffusivity of neighboring particles. Such studies might help 
to understand whether the effects shown here are consistent with an enhanced mobility 
(or reduced viscosity) of the surrounding fluid in the neighborhood of the tracer. We 

also plan to explore the possibility of optimizing tracer dynamics in deeply supercooled 
liquids, where the caging structure (and its dynamical consequences) are significantly more 
pronounced. 
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